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I treat the worldtube constraints which arise in the null-timelike initial-boundary value problem 
for the Bondi-Sachs formulation of Einstein's equations. Boundary data on a worldtube and initial 
data on an outgoing null hypersurface determine the exterior spacetime by integration along the 
outgoing null geodsics. The worldtube constraints are a set of conservation laws which impose 
conditions on the integration constants. I show how these constraints lead to a well-posed initial 
value problem governing the extrinsic curvature of the worldtube, whose components are related 
to the integration constants. Possible applications to gravitational waveform extraction and to the 
well-posedness of the null-timelike initial-boundary value problem are discussed. 
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I. INTRODUCTION 

It is extremely gratifying to contribute this article in appreciation of Josh Goldberg's friendship and guidance, 
5— I ' especially because this opportunity to recall Josh's early work has led me to an interesting approach to a current 
problem. Josh and I first overlapped in 1963 when he came to Syracuse University as a new Professor. At that time 
I was very busy finishing my PhD thesis and our interaction was fortuitous in more ways than one. Not only did we 
' share the same research interests but the next year, when I was looking for my first position. Josh recommended me to 
I> . the Aerospace Research Laboratory at Wright Patterson Air Force Base, where in his prior position he had organized 
CO ' a general relativity research group which included such budding young relativists as Roy Kerr. It was curious to 
me that the Air Force sponsored research in a topic with no apparent military relevance. Several years later, in a 
cost-cutting measure. Congress also came to this curious realization and the general relativity group was disbanded. 
However, during the intervening years the lab was a Camelot for basic research, much to the credit of Josh's legacy. 
It was there that I came to meet and work with the organizers of this volume, David Robinson and Ed Glass, two 
PhD students of Josh who came to the Lab on National Academy of Science postdoctoral fellowships. 

Much of Josh's research centered around conservation laws and null hypersurfaces [IHll- In the 1960's, these two 
topics were the focus of some of the most exciting results in gravitational theory. They came together in Bondi's |3| and 
Sachs's @ treatment of Einstein's equations via the characteristic initial value problem, which led to the formulation 
of the Bondi mass Mb and news function N as the quantities of prime physical importance for an isolated system. 
The conservation law 



((\N\^smed0d^, (1.1) 
du J 

which relates the retarded time derivative of the Bondi mass to the integral of the news function over the sphere at 
future null infinity, was the conclusive theoretical evidence that ended any serious debate over whether gravitational 
waves carry energy from a system. The conformal compactification of future null infinity into the boundary I"*" of 
an asymptotically flat spacetime put this mass loss relation into a well-defined geometrical setting Q. One of Josh's 
little known contributions to these results was the grant support of Bondi's King College group, which he arranged 
through the Lab. 

The Bondi mass loss equation is a result of certain constraints that arise in the characteristic formulation of 
Einstein's equations. The conservation law (jl.ip is obtained by applying these constraints at In this paper, I will 
discuss the content of these constraints when applied on a timelike worldtube of finite size but surrounding the matter 
sources of the gravitational field, as arises in the null-timelike initial-boundary value problem for the gravitational 
field. There is much overlap between my results presented here and Josh's work Q which supplies the basic ideas for 
interpreting these constraints as conservation laws. Josh applied these conservation laws to the theory of equations 
of motion for isolated systems being developed by Newman and his students 0, Q . Here I apply these conservation 
laws to a problem in numerical relativity that did not exist at that time. It is fitting, when I had a chance to talk to 
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Josh at a recent meeting on mathematical relativity, that he commented (essentially) "these are the same problems 
that we tried to solve in the 1960's" . That is true. The importance of those problems was recognized back then and 
only recently have many of them been elevated to a high level status by the mathematicians; and their importance 
to numerical simulation has also been recognized. But they are the same problems! 

In the Cauchy problem, initial data on a spacelike hypersurface Sq are extended to a solution in the domain of 
dependence I?(iSo) (which consists of those points whose past directed characteristics all intersect Sq). In the initial- 
boundary value problem (IB VP), data on a timelike boundary T transverse to Sq is used to further extend the solution 
to the domain of dependence I?(5o U T). 

The IBVP for Einstein's equations only recently received widespread attention due to its importance to numerical 
relativity foj, where the introduction of a finite artificial outer boundary is standard practice. It is essential for 
numerical evolution that the underlying analytic problem be well-posed, i.e. that the solution depend continuously 
on the data so that it is not destabilized by numerical error. The first well-posed IBVP was achieved for a formulation 
based upon a tetrad, connection and curvature as evolution fields [101] and subsequently for the harmonic formulation 
of Einstein's equations jTll |. However, much of the work in numerical relativity is based upon other formulations where 
the well-posedness of the IBVP remains an unresolved issue, incuding the nuU-timclike formulation. The properties 
of the worldtube constraints treated here is important for a clearer understanding of this characteristic IBVP. 

I begin with a short review of the nuU-timelike IBVP in Sec. [Til In the traditional approach to the Cauchy problem 
for Einstein's equations, initial data on a spacelike hypersurface with unit timelike normal are formulated in a 
purely 3-dimensional form in terms of the intrinsic metric /i^i, and extrinsic curvature k^^, of the initial Cauchy 
hypersurface. The components G^'^n^ of the Einstein tensor contain only first time-derivatives of the metric so that 
they constrain this initial data, i.e. the Hamiltonian and momentum constraints 

= (3)^^(fcM)2_^^^^M- (= 2G^,n''n'') (1.2) 

and 

= i^)V^{K-5^,kP) i^KG^pn"), (1.3) 

where '•'^-'Vp is the covariant derivative and '^■'i? is the curvature scalar associated with hf^^. Subject to these 
constraints, the Cauchy data determine a solution of Einstein's equations which is unique up to a diffeomorphism 

{d.m)- 

These constraints are elliptic partial differential equations which couple their solution with data on the boundary of 
the initial Cauchy hypersurface. Once these constraints have been satisfied initially, the dynamical equations used in a 
consistent hyperbolic reduction of Einstein's equations ensure that they continue to be satisfied in time in the domain 
of dependence of the Cauchy problem. Boundary conditions must then be formulated to preserve these constraints 
in the larger domain of dependence of the IBVP problem. That is the strategy which lies behind the constraint- 
free Cauchy evolution schemes used in numerical relativity. In characteristic versions of the IBVP, data are given 
on an initial null hypersurface u = const. This problem has the peculiar feature that the normal co-vector V^u is 
tangent to the hypersurface when interpreted as a vector k'^ — g'^^V^w by the standard technique of "raising indices" 
with the metric. As a result, the corresponding Hamiltonian and momentum constraints reduce to propagation laws 
consisting of ordinary differential equations (ODE's) along the characteristics of the null hypersurface. In the null- 
timelike version of the characteristic IBVP the worldtube integration constants for these ODE's uniquely determine 
the exterior spacetime. In SecHHwe review these ODE's and their integration constants. 



II. THE WORLDTUBE-NULLCONE PROBLEM 



A. The null cone formalism 



We use coordinates based upon a family of outgoing null hypersurfaces. We let u label these hypersurfaces, 
{A = 2,3) be labels for the null rays and r be a surface area coordinate. In the resulting x" = {u,r,x^) coordinates, 
the metric takes the Bondi-Sachs form 0, Q 

ds^ = - (e^'^— - r'^hABU^U^] du^ ~ 2e^^dudr - 2r^hABU^ dudx"^ + r^hABdx^dx^, (2.1) 



r 

where det{hAB) ~ det{qAB) — q, with qAB a unit sphere metric. The contravariant components are 

gvr ^ ^-2pY_ (2.2) 
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(2.3) 
(2.4) 
(2.5) 

where h-^^hsc = Sq. 

The Einstein equations G^jy — decompose into hypersurface equations, evolution equations and conservation 
laws. We express these equations following the formalism in The hypersurface equations correspond to the 

components G'j^Wi.u = and take the specific form 

= ^rh^''h^''hAB.rhcD,r (2.6) 
16 

(r^e^^^/i^Bf/,?),. = 2r^r-^(3^A)^^-r'h^''DchAB,r (2.7) 

2e-^^Vr = n- 2D^DAli - ^D^PDaP + r'^e'^^ DA{r''U^),r - ^r^e'*^ hAsU^^U^, , (2.8) 

where Da is the covariant derivative and TZ the curvature scalar of the 2-metric Hab- 

The evolution equations can be picked out by introducing a complex polarization dyad m'^ satisfying to^V^m ~ 
which points in the angular direction with components — (0,0, m"^) satisfying m^^fh^^ — h^^ , so that 
hABfTi^fh^ = 2, which determines up to the phase freedom m"^ — >■ e^^m^. The evolution equations correspond 
to the components m^^myG^i, = and take the form 

m^m''l^r{rhAB,u).r - ]^{rVhAB,r),r - 2e^^ {DaDbP + DaPDbP) 

+ hciADB){r^U^),r - \r^e-^^hAchBDU':,U^. + ^^r^hAB.DcU'^ 

+ r^U^DchAB.r - D^U^BhA)C..r + D^BU^hA)C,r^ = 0. 



(2.9) 



If we introduce the auxiliary variable Qa = r'^e HabU^ , the system of hypersurface and evolution equations can 
be cast into a hierarchy of first order ODE's, 

Ar- = J^pihcD) (2.10) 

{r^QA)^r = A/q(/icd,/3) (2.11) 

= Afu{hcD,l3,Qc) (2.12) 

Vr = MvihcD,f3,Qc,U^) (2.13) 

m^m''irhAB,u),r = Mh{hcD, P,Qc .U^ ,V), (2.14) 

where the A/'-terms on the right hand side can be calculated from the values of their arguments on a given 
u — const, null hypersurface. Each A/'-term only depends upon previous members of the hierarchy in the order 
{hcD, PtQctU''^ ,V). Thus, given Hab on the null hypersurface u = const., these equations can be integrated radi- 
ally, in sequential order, to determine /3, Qa, U^, V and m^m^ hAB,u on the hypersurface in terms of integration 
constants on an inner worldtube. For an inner worldtube given by r = R{u, x^), the necessary integration constants 
are 

/3U, Qa\b, U%., V\r, m^m^'hAB^ulR. (2.15) 
In addition, the location of the worldtube specified by R{u, x^) is another essential part of the data. 



III. THE BOUNDARY CONSTRAINTS AS CONSERVATION LAWS 



The components of Einstein's equations independent of the hypersurface and evolution equations are the conserva- 
tion conditions (called supplementary conditions by Bondi and Sachs) 

h^^'GAB = (3.1) 
G:, = (3.2) 
G: = 0. (3.3) 
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As was shown by Bondi and Sachs, the Bianchi identity 



V^Gi: = i^G^),^ + -g':^Gp. = (3.4) 



1 

-.V . 2. 

imphes that these equations need only be satisfied on a worldtube r = R{u^ x^). When the hypersurface and evolution 
equations are satisfied, the Bianchi identity for = r reduces to H^^Gab = so that p.ip becomes trivially satisfied. 
(Here it is necessary that the worldtube have nonvanishing expansion so that the areal radius r is a non-singular 
coordinate.) The Bianchi identity ioT: v = A then reduces to 

{r^G^A)^r = 0, (3.5) 

so that = if it is set to zero at r = i?(u, x^). When that is the case, the Bianchi identity for u = u then reduces 
to 

{r'Gl),r = 0, (3.6) 

so that GJ^ — also vanishes if it vanishes for r = R{u, x^). 

As a result, the conservation conditions can be replaced by the condition that the Einstein tensor satisfy 

eC^V.r - 0, (3.7) 

where ^'^ is any vector field tangent to the worldtube. This allows these conditions to be interpreted as flux conservation 
laws for the ^-momentum contained in the worldtube (is} . The unit norma to the worldtube A^^ lies in the direction 
V^(r-i?(u,x^)), i.e. 

N^^T^V^{r-Riu,x^)) (3.8) 

where ry is a normalization constant. Since we we are assuming that the hypersurface equations are satisfied, we can 
replace p.7|) by 

eO^N, = 0. (3.9) 



These are the boundary analogue of the momentum constraints (|1.3I1 for the Cauchy problem. In a treatment of a 
timelike boundary for the Cauchy problem, it has been pointed out [9] that the Cauchy momentum constraints (jl.3p 
must be enforced on the boundary. Here, in the characteristic initial-boundary value problem, it is the boundary- 
momentum constraints (|3.9p which must be enforced. 

Since ^^A^^ = 0, we can further replace p.9p by the condition on the Ricci tensor 

C^N^ = 0. (3.10) 

The Ricci identity 

^''i?;; = v^v^-^e"^ + v^vi'^e^i - vv^e" (3.11) 

then gives rise to the strict Komar conservation law 

V^yl"^^! = (3.12) 

when is a Killing vector corresponding to an exact symmetry. More generally, p. lip gives rise to the fiux conser- 
vation law 

P^{u2) - Pc(ui) = f ' dS^iW'W^e - V^V^"-^^)} (3.13) 



where 



P^^ j^dS^M"^^^ (3.14) 



and dS/^ii, and dSi, are, respectively, the appropriate surface and 3-volume elements on the worldtube. For the 
limiting case when i? — ^ cx), these flux conservation laws govern the energy-momentum, angular momentum and 
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supermomentum corresponding to the generators of the Bondi-Metzner-Sachs asymptotic symmetry group ^IS']. For 
an asymptotic time translation, they give rise to the Bondi mass loss relation 

Josh applied these conservation laws to a new treatment of equations of motions in general relativity Q. Here I 
pursue a different application to the mathematical basis of the worldtube constraints for the null-timelike IBVP. For 
this purpose it is useful to rewrite the conservation laws (jS.lOp in terms of the intrinsic metric and extrinsic curvature 
of the worldtube. 

The intrinsic metric of a worldtube embedded in the spacetime with unit spacelike normal iV^ is 

H^, = g^. - N^N, (3.15) 

and its extrinsic curvature is 

Kf,, = HP^pN,. (3.16) 

We then have the worldtube analogue of (|1.3p for the Cauchy problem, 

O^V^iKI^-Si^Kp i^H^G^pN"), (3.17) 

where 2?^ is the covariant derivative associated with i/^j/. These are equivalent to the conservation conditions (I3.10p 
and allow the conserved quantities to be expressed in terms of the extrinsic curvature of the boundary. For any vector 
field tangent to the worldtube, p.l7p implies 

2?p(r AT " eK) = ^^''^''^ (^A- - H^^uK). (3.18) 

In particular, this show that need only be a Killing vector for the 3-metric iJ^'^ to obtain a strict conservation law 
on the boundary. 



IV. THE WELL-POSEDENSS OF THE BOUNDARY CONSTRAINT PROBLEM 

The conservation conditions p.l7p constrain the boundary data for the nuUcone-worldtube problem. I now show 
that these constraints can be formulated in terms of a well-posed initial value problem intrinsic to the worldtube. 
The traditional 3-1-1 decomposition of spacetime used in the Cauchy formalism is not applicable to the nuUcone- 
worldtube problem since a foliation by hypersurfaces with a null normal gives rise to a degenerate 3-metric. However, 
an analogous 2-1-1 decomposition can be made on a timelike worldtube. Let vector and tensor fields intrinsic to the 
worldtube be denoted by w", etc. 

Let be an evolution vector field on the worldtube, i.e. the flow of carries an initial spacelike cross section t = Q 
into a t-foliation of the worldtube, with Ct^t = 1. Coordinates of the initial slice then induce adapted coordinates 

= {t,y"^) on the worldtube by requiring Ct<^y^ — 0. Thus the choice of i° and initial coordinates y^ fix the gauge 
freedom on the worldtube. 

The intrinsic 3-metric Hab of the worldtube has the further 2-^1 decomposition 

Hab = -TaTb + R^hab, H"^ = -T^T^ + R'"^ h"^ . (4.1) 

where is the future timelike unit normal to the foliation and R^hab is the intrinsic 2-metric of the t — const 
slices, with h^^hsc = ^c- Here, for later application to the characteristic problem. We have introduced the surface 
area factor R{t,y'^), so that det{hAB) = det((5^B) = Q, where Qab is a unit sphere metric. For convenience, we 
chose stereographic coordinates y^ = {q,p) for the unit sphere metric for which QAsdy^dy^ = Q^^^{dq^ + dp^) with 
Qi/2 —4/(1-1- (jf2 +p2-j2_ gjjjQg j^AB j^g^g (-|-4_) signature, it can be put in the matrix form 

where 7 and a represent the two degrees of freedom. A specific choice of polarization dyad associated with this 
representation is 

— Q~^^^ ^e~^(coshQ! — i sinh a), ie'''(cosha -|- i sinha)^ . (4-3) 

The operator 



hi = St + T^Tb 



(4.4) 
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projects tensor fields into the tangent space ortliogonal to T°. The relations 

= AT" + 3" , Ta^ -Adat (4.5) 

are the 2 + 1 relations defining the lapse A and shift i?" on the worldtube analogous to the standard 3+1 decomposition 
of the Cauchy problem. In the adapted coordinates, Hab — R^hAB is the metric of the 2-surfaces of constant t on 
the worldtube, 

A^ = -Htt + HabB-^B^ (4.6) 

is the square of the lapse function and 

= R-^h^^HtB (4.7) 
is the shift vector. We have Ta — —Adat, the contravariant components of the worldtube metric are 

jjtt ^ _^-2 

jjtA ^ A-^B^ (4.8) 
H^"" = R-^h^^ -A-^B^B"". 

and we again have the dyad decomposition h"^^ — m!^°'rffi with m" = (0,m'^). 

The mathematical structure of the conservation conditions (j3.17p is simplest to analyze in a worldtube gauge in 
which the lapse A = 1 and the shift B^ — 0. This corresponds to the introduction of Gaussian normal coordinates 
on the worldtube, which is always possible locally. In this Gaussian gauge, Hu = -ff** = —1, HtA — H*^ — 0, 

Let Va denote the 3-connection on the worldtube associated with Hat- Then in the Gaussian gauge, the conservation 
conditions p. 171) reduce to 

V.iKl - 5lK) = -j^d, [r^VQ{kI - S'aK)^ + ii?,?^(ifBC - HbcK) = (4.9) 

with components 

-^dtih''^'KBc)-VBK^ -^H^^{Kbc-HbcK) - (4.10) 

-^dt{R^KAt)-VBK^+VAK = 0. (4.11) 

These equations can be re-expressed in terms of the 2-connection cIa on the slices of the worldtube associated with 
hAB as 

dt{h^''KcD)-dA{h^''KtB) = lh^t''iKBC-hBch''''KDE)~^R.t{h''^KBc-2R^K) (4.12) 

dtih^^KtB) - In^^'dBih^'^KcD) = R'H^f^KtB + h^^dBiK^ - Ih^KE) - h^^dBK 

+ ^R,B{h''''K^-h^''KS). (4.13) 

We interpret (|4.12p - (I4.13P as a system of equations for h^'~' Kbc and Kai given their right hand sides. As such, 
they take the form of a symmetric hyperbolic system (cf. flG^). This requires that, up to lower differential terms, the 
left hand sides have the form 

dtV - M^^^dAVf^ (4.14) 

where, for each ^-component, T-L^^Mp^ is symmetric in (7,/?) for some symmetric, positive-definite matrix Tija- In 
the present case, we have = ^(y*,F«,yP) = {h^^ KcD,h'i^ Ket^hP^ KBt)- The component is given (up to 
lower order terms) by 



2 

qqq 

mp 



M"'^ = i ( I . (4.15) 
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Referring to (|4.2[) . the symmetrizer is 

/ (2i?2cosh2a)-i \ 

n-ya^i e^T tanh2Q . (4.16) 

\ tanh2a e'^'^ J 

The component M^^ has the same symmetrizer. 

Such symmetric hyperbohc systems have a well-posed Cauchy problem, i.e. there exists a unique solution which 
depends continuously on the values of h^'~^ Kbc and KtA on the initial slice of the worldtube. In addition to lower 
order terms in y", the right hand sides of (I4.12p and (|4.13l) depend upon Kg - {l/2)h^K'§, K, Hab and R. Here 
Kg — {l/2)hgK^ is determined by ■m°-m!^Kab or, in tensorial form, by — ^hlh'^)K^. Thus, in the Gaussian 

gauge, we have established the following Worldtube Theorem: 

Given Hab, mf^m^Kab and K , the worldtube constraints constitute a well-posed initial-value problem which determines 
the remaining components of the extrinsic curvature Kab- 

This theorem extends to any gauge. More generally, the conservation conditions p.l7p take the form 

T^db{hiK^)~dB{h^K^T'') = St (4.17) 
T^db{h\KiTd) - \dA{hiK^a) = Sa, (4.18) 

where the source terms St and Sa are determined by Hab , 'm°"m^Kab, K and lower order terms. The system (|4.17p - 
(|4.18p can again be symmetrized. 

It is important to note that whether or not Hab, mf^m^Kab and K can be chosen independently depends upon the 
choice of gauge conditions and formulation of Einstein's equations. In the next Section, we discuss this issue in the 
context of the Bondi-Sachs formulation. 



V. CONSERVATION LAWS AND THE CHARACTERISTIC INTEGRATION CONSTANTS 

The worldtube theorem established in Sec. IIVI formulates the boundary constraints as a well-posed initial value 
problem given the required source terms. We now explore the implications for the integration constants necessary to 
integrate the characteristic hypersurface and evolution equations. The integration constants (|2.15l) constitute 8 real 
functions on the worldtube. The areal radius of the worldtube R{u,x^) is a 9th function necessary to determine a 
unique solution. Thus no more than 9 pieces of worldtube data can be freely specified. This matches the number 
of functions assumed in the worldtube theorem, i.e. -ffab, mf^m^Kab and K . However, the boundary constraints 
introduce 3 relations between the 9 pieces of worldtube data. 

We denote by y° — (t, y^) the coordinates intrinsic to a foliation of the worldtube by topological spheres t = const, 
with angular coordinates y^. The propagation of these coordinates along the null geodesies of the outgoing u = const 
null hypersurfaces, which are uniquely determined by this foliation, induces Bondi-Sachs coordinates x^^ = (u, r, x^) in 
the exterior spacetime with a;" = {u,x^) — on the worldtube. When convenient we will switch from 4-dimensional 
notation with coordinates to 3-dimensional notation with coordinates y". We have 

^ — ^ + i? ^ (5 1) 

where we continue to denote partial derivatives by commas when it does not lead to ambiguity. In addition, in this 
Section, we denote d/dy"" = da- 

The metric on the worldtube r = R[u,x^) which is induced by the Bondi-Sachs metric ()2.1|) is 

Habdy'^dy*' = -e^/^^dt^ - 2e^'^ R^adtdy"" + R^hABidy"^ - U^dt){dy^ - U^dt). (5.2) 
R 

Here V{t,y^) is related to the Bondi-Sachs variable V{u,r,x^) by V{t,y^) = V{t,R{t,y^),y'^). Similarly, 
hAB{t,y^) — hAB{t, R{t,y^),y^), etc. We again denote by dA the 2-connection on the slices of the worldtube 
associated with Hab- It is important to distinguish between dA and the 2-dimensional covariant derivative Da as- 
sociated with hAB on the (r = const, u = const) Bondi-Sachs spheres. In the tangent space of the worldtube, we 
again have the dyad decomposition /i"'' = m'^°'fh!'^ with m" = (0,m'^). It is also important to distinguish between 
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vectors tangent to the worldtube and vectors tangent to the r = const Bondi-Sachs coordinate surfaces. In the 
Bondi-Sachs coordinates , a vector tangent to the worldtube has components = {v"',v^) where = V^R.a- 
Thus mf'df, = m'^dA + m^R,Adr. This same procedure aUows us to express and in terms of Bondi-Sachs 
components and iJ'"^, e.g T^'S^ — T'^da + T'^R^adr- In this regard, it is useful to recall that forms pull back and 
vectors (and other contravariant fields) push forward under the embedding map from the worldtube to the spacetime. 

In order to interpret the implication of the boundary constraints for the nuUcone- worldtube problem we now relate 
the intrinsic metric and extrinsic curvature of the worldtube to the integration constants (j2.15|) for the hypersurface 
and evolution equations. From (|4.6p and (14. 7p . the lapse and shift corresponding to (|5.2p are given in terms of 
Bondi-Sachs variables by 

A2 = e^^i^ + 2Rt) + R^hAB{B^ + U^){B^ - U^) (5.3) 



and 



= -R-^e^^h^''R,B - U^- (5.4) 



The geometric properties associated with the t-foliation of the boundary, with normal Ta = —Adat, are the 3- 
acceleration 

^ T'^VaT'' (5.5) 
and the 2-dimensional extrinsic curvature of the i-foliation 

Kab = KV.Tt. (5.6) 

The independent components of and Kab are 

m^ab = m^dB log A (5.7) 

/i^^^afc = -A-' (^dciR^B^) ~ dt{R^)^ (5.8) 

and 

R^ R^ 

m'^m'Kab = -z-rm^m^dthAB -r-m^niDdcB^ . (5.9) 

The 3-dimensional extrinsic curvature of the worldtube is K^^ ~ H^Va^iy where the unit outward normal is given 

by 

N^ = r,'^^,{r-R). (5.10) 
The normalization condition expressed in Bondi-Sachs coordinates gives 

7/-' - -I- 2e-^f'R,u + 2e-^"U^R.A + R-^h^"" R,aRb = e'^^ A\ (5.11) 
R 

so that -q — e^^ A^^ . The outgoing null vector normal to the i- foliation of the worldtube is 

In Bondi-Sachs coordinates, K^dfj, is proportional to dr- The proportionality constant is determined by the normal- 
ization condition K^N^ = 1, which evaluated on the worldtube gives 

K^'d^ = r]-^dr = e-^I^Adr. (5.12) 
The computation of i^^j^ in terms of Bondi-Sachs coordinates can be simplified by the construction 

= H^V^N, = H^^H^^V^Np (5.13) 

= H^X.)^'^Kp - H^H^^^V^Tp (5.14) 

= - /i^/i^23„r/3 + /i^^r,)r^2?^r„ (5.15) 
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so that 

1 



Kf,r, = -H^HPCK9ai3 - K^i. + a(^T^). (5.16) 

Expressions for the components of CKgap in terms of the Bond-Sachs variables, which are vahd for an arbitrary 
gauge, are given in Appendix]^ Collecting the pieces of (|5.16p . the components of the extrinsic curvature are 

h^^Kab = 2Ai?e-2/3 + ^-1 (dc{R^B^) ~ dt{R^)] (5.17) 



and 



1 r>2 7-)2 

m''m}'Kab = -AR^e-'^^m'^m^hAB,r - Trrm'^m^ dthAB + -r-rrf modcB'^ (5.18) 

m'^T^'Kab = -Im^QA - m^R^.r " 7^) + \m^m''hAB,rrff R,c - dc \og{e-^^ A^) (5.19) 
Z ZH 4 Z 



- + -^^^h^''R,AR.B - T'^da logie-^f'A). (5.20) 



The interplay between the boundary constraints on Kab and the characteristic integration constants is complicated 
by the choice of gauge and the choice of free data. We consider two complementary scenarios which simplify the 
discussion of the underlying problems. 



A. Waveform extraction 



In the first scenario, the Bondi-Sachs integration constants {f3,V,U^ ,QA,hAB,R) are obtained from the metric 
in the neighborhood of the worldtube which is supplied by the numerical results of a 3 + 1 Cauchy evolution. This 
provides the inner boundary data for a numerical characteristic evolution on a Penrose compactified grid, which yields 
the waveform at T+ . This approach, called Cauchy-characteristic extraction [13, [13] , is used in numerical relativity to 
obtain the waveform at 1+ without the near field ambiguities introduced by a finite outer boundary for the Cauchy 
evolution. In this application, error in the Bondi-Sachs integration constants results from numerical error inherent in 
the Cauchy evolution, from error in computing the Jacobian from the 3 4-1 Cartesian coordinates to the Bond-Sachs 
coordinates and error from interpolation onto the worldtube. This feeds into a corresponding error violation of the 
worldtube conservation laws, i.e. constraint violation. 

Alternatively, the Cauchy data can be used to provide {Hab, Tn°''m^Kab, K) on the boundary. The worldtube theorem 
can then be applied to determine the remaining components of the extrinsic curvature {h^K^,m'^TbK^) using the 
well-posed evolution system provided by the conservation laws. This approach could be used to control the constraint 
violation in the characteristic evolution which results from error in the Cauchy data. 

The implementation would proceeds as follows. Since {A,B'~' ,R) are provided from Hab, the component h'^K^, as 
given in (|5.17p . determines the integration constant for (3. The shift B'-' then determines the integration constant for 
via (|5.4p : and the lapse A then determines the integration constant for V via (15. 3p . Next, wf'TbK'^ determines 
the integration constants for in the following way. From m'^m^Kab, as given in (|5.18p . we can determine a := 
jm^^m^drhAB, the optical shear of the null hypersurfaces emanating outward from the worldtube. Given ct, we 
can determine dr(3 using the hypersurface equation (|2.6p . The integration constant for is then determined from 
mf'TbK''^, as given by (lETO)) . 

We have thus shown that the worldtube theorem leads to the following Extraction Corollary: 

Given Hab, m'^m^Kab and K, the worldtube constraints constitute a well-posed initial-value problem which determine 
the the Bondi-Sachs integration constants {13,V,U^ ,QA:hAB: R)- 
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B. The initial-boundary value problem 

In the forgoing waveform extraction scenario, the nine integration constants {f3,V,U^,QA,hABTR), which are 
required to integrate the Bondi-Sachs equations were supphed by a Cauchy evohition in a manner consistent with 
Einstein's equations. The data consisted of the boundary metric Hat, (6 functions), m°'m^Ka,b (2 functions) and K 
(1 function). The constraints were then enforced via the worldtube theorem to determine the remaining components 
of Kab^ which in turn supphed the integration constants. 

In an initial-boundary value problem, boundary data consistent with the constraints must be prescribed apriori, 
i.e. before the evolution is carried out. Enforcement of the boundary constraints has been the major difficulty in 
attempting to show that the various formulations of the gravitational initial-boundary value problem are well-posed. 
In this regard, the characteristic formulation is no exception. The coupling between the Bondi-Sachs evolution system 
and the boundary constraint system is complicated. The details depend upon the choice of free boundary data and 
the choice of gauge conditions adopted on the boundary. We illustrate this problem with two examples. 

1. Constant R boundary data 

Consider first the case in which the the 5 worldtube integration constants for {P,U^,hAB) are prescribed freely 
along with a constant value of the areal radius R. We then attempt to prescribe the remaining 3 integration constants 
for {V,Qa) via the 3 boundary constraints (I4.17|) and (j4.18p . For simplicity, assume the boundary data are /3 = 
U"^ = 0, hAB = Qab (where Qab is the unit sphere metric), along with R — const. The lapse (j5.3p and shift (|5.4|) 
corresponding to this data reduce to 

A^^V/R, B^=0. (5.21) 

Along with the initial data at i = 0, 

/iAB(0,r,a;^) = Qab, V{Q,R,x^) = R-2M , Qa{0, R, x^) = 0, (5.22) 

this worldtube data determine a mass M Schwarzschild spacetime in spherically symmetric Bondi coordinates. 

Now, with the same boundary data, let the initial data hAB{0,r, x^) consist of a pulse whose support is isolated 
from the worldtube. This evolves to produce ingoing radiation so that the spacetime in the neighborhood of the 
boundary is no longer Schwarzschild in the domain of dependence of the initial pulse. The boundary constraints for 
this problem reduce to 



dtK§~-dB{AKf)^0 (5.23) 



1 
A 

and 

dtijKBt) - ^dBiAKS) = dc(^A{K'^ - ^S'^kE)^ - dB{AK) + kScIbA, (5.24) 

where the components of the extrinsic curvature are 

h^^KBc = '2AR (5.25) 

m^T^'Kcb = -m^i^Qc + jdcA) (5.26) 
AR"^ 

m^m^KBc ^ -^m^rn^hBC.T (5.27) 

and 

'lA 1 1 

Before the pulse hits the worldtube, V{u, R, a;^) = i? — 2AI and A^ = 1 — 2M/R, but afterward the worldtube values 
of hscr and Qa become time dependent. 
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If these boundary constraints were to constitute a well-posed problem for Kg and Kj^ then they would supply the 
boundary values of A and Qa necessary to determine the remaining Bondi-Sachs integration constants. However, two 
problems arise from the right hand side of (j5.24|) . First, the term — \S'^K^ is determined by the optical shear 
a at the boundary. But this shear is not known until the evolution is carried out and it depends upon the detailed 
shape of the initial radiation pulse. Thus the boundary system (15.231) - (|5.24l) is coupled to the evolution system. 

The second problem is more serious. For this system the evaluation of the hypersurface equations on the boundary 
implies 

P'- = -^K''^AB,r^^(J& (5.29) 

and 

Vr = l + \h^°dcQD - \h^°Qc,rQD,r (5.30) 

3/4/47? 11 1 1 

^ = 2i? + + 2Ai?(^ + ^h^^'dcQD ~ -^h^^Qc.rQo.r) + ^dtA. (5.31) 

As a result, in addition to coupling the boundary and evolution systems, the term dsiAK) on the right hand side of 
()5.24|) introduces a h'-^^dsdcQD term which alters the principle part of the boundary system so that it is no longer 
guaranteed to determine A and Qa in a well-posed manner. 



so that 



2. Trace K boundary data 

The problem in the preceding example arises because the boundary data for R and K (the trace of the extrinsic 
curvature of the boundary) cannot in general be given independently, as is required for application of the worldtube 
theorem. Specification of i? = const boundary data simplifies the boundary system since (15.11) then reduces to 
d/dy"" — d/dx"-. However, this simplification is offset by the complicated way in which K changes the principle part 
of the boundary system and affects the well-posedness through the dsiAK) term in (|5.24p . 

The only way this complication can be avoided is to prescribe K{t^ y^) as explicit boundary data, so that it does 
not affect the principle part of the boundary system. In Bondi-Sachs coordinates, R and K are related by 

K = Ht^'^V^N, = r^m^d^d.ir - R) - Ht'TP^^Np = -rjH^'>dadi,R - iJ^^^F^^iVp. (5.32) 

Here the timelike nature of the worldtube ensures that H°-''dadb is a wave operator. Thus, given K{t,y^), the lapse 
A, the shift B^, the conformal 2-metric h"-^ and the Bondi-Sachs Christoffel symbols H'^'^T'^^Np on the worldtube, R 
can be determined from its initial data by a well-posed quasilinear wave problem based upon (15 32p . This relationship 
between the function locating the worldtube and the extrinsic curvature scalar of the worldtube was first pointed out 
in the Friedrich-Nagy [lo| treatment of the initial-boundary value problem. 

Although this might at first sound like a promising approach, it leads to serious difficulties. Foremost, if R^a 
then the boundary constraint (|4.18p couples the evolution system with the boundary system in a way which makes 
the formulation of a well-posed evolution-boundary system appear to be intractable. This arises from the terms 

-■m^R^Al3,r + ^m^m^ hABj-m^ R.c 

in the expression (|5.19p for mf^T^Kab- Here and hAB,r cannot be determined without knowledge of the evolution. 
Consequently, the boundary constraint governing m°'T^Kab cannot be used to determine the Bondi-Sachs integration 
constant for Qa independently of the evolution. 

One possible way to circumvent this problem would be to pick a gauge for the boundary in which R a = 0, i.e. 
R — R{t). Since i? is a scalar density defined by the determinant of the 2-metric of the boundary slices, this can be 
achieved via the Jacobian of an appropriately chosen angular transformation. However, this raises the new problem 
of how to pick explicit data for K which would be consistent with R_a = 0. 



VI. DISCUSSION 



We have shown how the boundary constraints for the Bondi-Sachs equations can be posed as a symmetric hyperbolic 
system governing the evolution of certain components of the extrinsic curvature of the worldtube, as described by the 
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worldtube theorem in Sec. IIVI In Sec.|V]we described how these extrinsic curvature components were related to the 
integration constants for the Bondi-Sachs system. The apphcation of the worldtube theorem requires knowledge of 
the intrinsic metric of the worldtube and the remaining components of the extrinsic curvature. We considered two 
different versions. 

The first application was to waveform extraction. In that case, the data {Hab, m°'m^Kab, K) necessary to apply the 
worldtube theorem are supplied by the numerical results of a 3 + 1 Cauchy evolution. The remaining components of 
the extrinsic curvature can then be determined by means of a well-posed initial value problem on the boundary. The 
integration constants {13,V,U^ ,Qa, hAB, R), for the Bondi-Sachs equations are then detemined. This approach can 
be used to enforce the constraints in the numerical computation of waveforms at 1+ by means of Cauchy-characteristic 
extraction. 

In the second application, we considered the initial-boundary value problem, for which boundary data consistent 
with the constraints must be prescribed apriori, i.e. independent of the evolution. The object was to obtain a well- 
posed version of the characteristic initial-boundary value problem. However, the complicated coupling between the 
Bondi-Sachs evolution system and the boundary constraint system prevented any definitive results. Two choices of 
free boundary data and boundary gauge conditions were explored. In both cases, the Bondi-Sachs choice of areal 
coordinate r complicated the analysis. This results from the way that the coordinates and geometry are mixed, i.e. 
r cannot be assigned freely on the boundary without specifying its area. It is possible that other formulations of 
the characteristic initial-boundary value problem might be more amenable. Bartnik 19i] previously explored a quasi- 
spherical version, in which the 2-metric Hab is transformed into a conformally unit sphere form. He found similar 
complications in trying to establish well-posedness as in the Bond-Sachs case. A formal computational algorithm 
for the evolution-constraint system was possible, but the well-posedness of the corresponding initial-boundary value 
problem, which is necessary for numerical stability, was not clear. 

Another possibility is to choose a gauge in which the areal coordinate r is replaced by an affine parameter A, so that 
the affine freedom allows the specification A = on the boundary, independently of its geometry. Rendall [1^ has 
shown that such a characteristic initial-boundary value is well-posed in the double null case where the boundary is 
also a null hypersurface. However, Kendall's approach cannot be applied to the corresponding null-timelike problem. 
Although the full treatment of the null-timelike problem lies outside the scope of this paper, it is clear that the 
wordltube conservations laws must enter in an essential way. 

Appendix A 

We use dug), 

K^'d^ = e-^^Adr, (Al) 
to simplify the calculation of CKgai3 in Bondi coordinates. We have 

CkQcP = e-^PAgo^p^r + e^P [K^dpie-^l" A) + KMe-^^A)^ . (A2) 

Thus 

H^H^£K9ap = e-2^Ai?^i?,^5a/3,r + c^^ A-^ (^K^Op {e'^^ A) + Kpd^{e-^^A) 

= e-^^AH^HUg^p^r + T^V, logie'^f^A) + T,V,, log(e-2^*A). (A3) 



In Bondi-Sachs coordinates, the tangents to the worldtube foliation have components = (0, to^/J^b, Jti"^) and 
T^' = A-i(l, R^t - B'^Rc, -B^). This leads to the components 

m''m''H^H[;CKgap = -iARe-^^ (A4) 
m>'T7fH'^H^CKgap = AR'^e-^'^m^m^hAB,r (A5) 

m^T''H^H!^CK9c.(S = -tii^'Qa - 2M^R^AP,r - e-'^^m'' qabAU^ + B^) - m^d^ \og{e~^^ A) 

= -m^QA - 2m^R,A{f3,r - 4) + \m'^m^hAB,rfn^ Rc - m^dA log(e"2;3^) ^^g^ 

R 2 
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and 




- -±^e'0hfR,AR,B + -^e'^h^''R,AR,B - 2T«a„ log{e-'^A). 



(A7) 
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